Coexistence of Single and Double-Quantum Vortex Lines 
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We discuss the configurations in which singly and doubly quantized vortex lines may coexist in 
a rotating superfluid. General principles of energy minimization lead to the conclusion that in 
equilibrium the two vortex species segregate within a cylindrical vortex cluster in two coaxial domains 
where the singly quantized lines are in the outer annular region. This is confirmed with simulation 
calculations on discrete vortex lines. Experimentally the coexistence can be studied in rotating 
superfluid 3 He-A. With cw NMR techniques we find the radial distribution of the two vortex species 
to depend on how the cluster is prepared: (i) By cooling through T c in rotation, coexistence in the 
minimum energy configuration is confirmed, (ii) A glassy agglomerate is formed if one starts with 
an equilibrium cluster of single-quantum vortex lines and adds to it sequentially double-quantum 
lines, by increasing the rotation velocity in the superfluid state. This proves that the energy barriers, 
which separate different cluster configurations, are too high for metastabilities to anneal. 



PACS: 67.57.Fg, 52.25.Wz, 74.60.Ge 



I. INTRODUCTION 

Superfluid 3 He-A is the only known quantum sys- 
tem where topologically stable vortex lines with different 
quantization appear simultaneously. The doubly quan- 
tized vortex (DQV) is formed at low critical velocity and 
is usually obtained when superfluid 3 He-A is accelerated 
to motion, eg. in a rotating container jjj. The singly 
quantized vortex (SQV), in turn, may have lower energy, 
depending on the magnitude of the magnetic field and 
flow velocity, but because of its large critical velocity it 
is in practice only formed on cooling through the super- 
fluid transition in nonzero flow §. 

Normally one will find either double or single-quantum 
vortex lines in a rotating container, depending on how the 
state has been prepared. In the homogeneous rotating su- 
perfluid, transitions between different vortex textures are 
of first order and therefore in principle sharp. However, 
topological stability forbids a transition from one type of 
existing vortex lines to another. Thus SQV and DQV 
lines, once they have been formed, will maintain their 
quantization in the rotating container, as long as they 
are not removed by deceleration and annihilation on the 
container wall or by warming up to the normal phase. 
This metastability can be exploited to prepare a rotating 
state where both types of vortex lines are simultaneously 
present. 

A priori it is not clear what configuration a rotating 
state with two different vortex species would take: Do 
the arrays always have the same configuration or does 
this depend on the process by which the state has been 
created? Similar questions have been studied in the case 
of point-like objects, namely trapped ions in a 2 or 3- 
dimensional point charge system ||. We investigate the 
coexistence of SQV and DQV lines both theoretically and 



experimentally. We study the free energies of different 
vortex-array configurations, which consist of a mixture of 
the two species or of segregated domains. We apply both 
the continuum model and use numerical simulation of dis- 
crete vortex lines to show that a configuration with SQV 
vortices surrounding an inner central cluster of DQV is 
lower in energy than any other configuration. The en- 
ergy difference is small and disappears in the continuum 
approximation. 

Axially symmetric configurations of coexisting SQV 
and DQV lines are observed in our measurements when 
liquid 3 He is slowly cooled at constant rotation from the 
normal phase to 3 He-A. Our observations of these states 
agree with the theoretical predictions: the SQV lines are 
found in an outer annulus which surrounds a central do- 
main of DQV lines. Segregated configurations, which do 
not have full axial rotation symmetry, can also be pre- 
pared. This is accomplished by adding DQV lines to 
an existing equilibrium SQV cluster, by slowly increas- 
ing the rotation velocity at constant temperature in the 
superfluid state. The resulting configuration does not 
correspond to the global energy minimum, but rather re- 
sembles a glassy agglomerate. 



II. ROTATING HE-A. 

Rotating vortex states are driven and stabilized by su- 
perfluid counterflow (CF) , the difference in the velocities 
of the normal and superfluid fractions, v = v„ — v s . The 
equilibrium states are found by minimizing the free en- 
ergy functional j| , which controls the spatial distribution 
of the order parameter. On the macroscopic scale, vor- 
tex textures in 3 He-A are rectilinear structures, oriented 
parallel to the rotation axis, which extend from the top 
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to the bottom of the rotating container. These vortex 
cells are confined by the Magnus force to a regular array, 
a cluster of vortex lines. The cluster is coaxial with the 
cylinder || and isolated from the cylindrical wall by an 
annular layer of vortex- free CF. Within the cluster the 
areal density of linear vortex cells has the solid-body- 
rotation value n — 20,/ k. Here k = vkq is the circulation 
of a vortex line with the quantum number v and the cir- 
culation quantum kq = h/(2m^) = 0.066 mm 2 /s. 

In a vortex cluster of infinite diameter a transverse 
cross section would display a perfect 2-dimensional lat- 
tice with triangular nearest-neighbor coordination, which 
arises from the inter-vortex repulsion. In a finite cluster 
the central region is expected to display a triangular lat- 
tice, which approaches crystalline order, but towards the 
outer boundary it is deformed by the surrounding annular 
CF more and more towards a configuration of concentric 
rings of vortex lines [|l0| . Due to the competition between 
volume and surface interactions long range order is thus 
relatively poor in finite-size arrays at the usual exper- 
imentally accessible rotation velocities. Here the total 
number of lines is typically a few hundred, the radius R 
of the container a few mm, the radius of the Wigner- 
Seitz (WZ) cell of the vortex lattice ro ~ 100 /xm, and 
the vortex core radius £ ~ 10 pan fTl]] . 

The order-parameter texture within a unit cell of the 
periodic vortex lattice has been derived from numerical 
minimization of the free energy as a function of the ro- 
tation velocity CI and the applied field H f^M[|. These 
calculations reveal a complex phase diagram @||, with 
quantized vorticity organized in several topologically and 
structurally different vortex textures jTj. SQV vortex 
textures, which include a singular vortex core, have only 
been described approximately, by solving for the order 
parameter distribution with analytical trial functions . 
Of these different vortex textures we shall here only be 
concerned with two, the continuous DQV and the singu- 
lar SQV. 

In an array consisting of DQV lines the inter-vortex 
distance is larger by \[2 than in one formed from SQV 
lines. With coexisting SQV and DQV lines this large dif- 
ference opens the question whether sufficient mismatch 
and distortion is created in the nearest-neighbor coordi- 
nation that a disordered glassy state results. Or could 
it perhaps lead to a "chemical" structure where both 
species occupy fixed positions within each unit cell of 
a periodic vortex lattice. A third possibility is an ar- 
ray consisting of segregated domains of one species only. 
These are the questions which we shall discuss next. 



III. ENERGY MINIMIZATION 



A. Continuum approximation 

On the macroscopic level the energy difference between 
SQV and DQV arrays depends on the difference in the 
quantization numbers and inter-vortex distances, if the 
latter are much larger that the core sizes. Here we ne- 
glect the smaller energy contributions associated with the 
structure of the vortex core. Consider the free energy of 
the rotating superfluid 
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(v 2 - 2v s v rot )d 2 r , (1) 
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where v rot = f2 [z x r] , is the rotation velocity, and 
L the angular momentum of the fluid. We express the 
superflow velocity v s in the continuum approximation 
and assume the logarithmic approximation for the vortex 
line energy. In the laboratory frame we then have |12| 
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The second term in the integral represents the energy 
of a vortex line. If different vortex species are present, 
the circulation K a — v a KQ, the inter- vortex distance r a , 
and the core radius £ a become functions of the spatial 
coordinates (here a = 1 or 2). 

Let us consider a configuration with N± SQV lines and 
N2 DQV lines. The overall radius R v of the vortex cluster 
is determined by the total circulation: 
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In the stationary state the areal densities of vortex lines 
are proportional to their circulation, 
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The average velocity field is v s = v s0 where 
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if r > R v . 
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The free energy in Eq. (|^) is now given by 
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where E = np s n 2 R 4 , S = ttR 2 , and S a = K a N a /(2n) is 
the area occupied by SQV or DQV lines, such that Si + 
52 = ttR 2 ,- The normalized vortex- line energy densities 
in the logarithmic approximation (see Eq. (p9|)) arc 
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These are proportional to the first power of distinct 
from the single- vortex energy which is proportional to k 2 
(where a = 1 or 2). 

These expressions are valid when r a 3> £ a HI- The 
energy densities of SQV and DQV arrays become equal 
when 
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The rotation velocity for the transition becomes 

2k q £ 



n* = 
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To explore the consequences from the continuum ap- 
proximation further, we minimize the energy in Eq. (|^) 
with respect to variations of both Ni and N2 . This gives 
the equilibrium configuration, which one would expect 
to find after an adiabatic transition at constant rotation 
from the normal to the superfluid state. Several conclu- 
sions follow: 

(1) If £2 > £i, only SQV lines remain since S2 shrinks to 
zero (Eq. (0)). In the opposite case, Si shrinks to zero, 
and only DQV lines remain. The DQV becomes more 
favorable, when f2 > f2*. According to our measurements 



this occurs above 0.6 rad/s (Sec. VAl). 

(2) The radius R v of the cluster has its equilibrium 
value when the width of the surrounding vortex-free CF 
annulus becomes 



^a(eq) 



= R^7j2. 



(9) 



where the outer boundary of the cluster consists of type 
a vortices. 

(3) The energy in Eq. (^|) does not depend on the con- 
figuration of the vortex cluster with given Ni and N2- 
Consequently, the continuum approximation in the loga- 
rithmic limit does not discriminate between different vor- 
tex cluster configurations. 



B. Discrete vortex lines 

1. Sequentially formed arrays 

The continuum approximation in Eq. (|J) does not in- 
clude the surface energies at the interfaces between the 
vortex cluster and the vortex-free region or between seg- 
regated domains occupied by different vortex species. 
Nor does it account for interactions between individ- 
ual vortex lines and the container wall. To identify 
configuration-dependent contributions to the free energy, 
we need to look at vortex clusters composed of discrete 
vortex lines. We perform numerical calculations on vor- 
tex arrays which consist of Nj discrete rectilinear vortex 
lines with coordinates Tj, circulations Kj = VjHo, where 



1,2, and core radii The superfluid velocity 



produced by the lines is 
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The last term accounts for the image vortices with coor- 
dinates r'j — TjR 2 /^, of which each one corresponds to 
one true vortex located at rj. 

The interaction potential for vortex lines in the rotat- 
ing superfluid is equivalent to the Coulomb potential be- 
tween charged lines immersed into a uniform background 
charge of opposite sign. The density of the background 
charge is proportional to Q. The potential acting on each 
line is the sum of the potentials of the background charge, 
of the direct logarithmic Coulomb interline interaction, 
and the interactions with the wall reflections. 

The equation describing the dynamics of vortex lines 
has the form 



= (l-b')V s -b[. 
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where Yj is the velocity of the jth vortex line and v' s the 
superfluid velocity produced by all other vortices at the 
position of the jth reference line, while b and b' are the 
reduced Hall-Vinen mutual friction parameters Q3j. 

To illustrate the dynamic simulation in vortex array 
formation, we consider a standard experiment when a 
vortex array is formed in the superfluid state by increas- 
ing slowly the rotation velocity linearly from zero Q. In 
this process one vortex line is formed after another pe- 
riodically, every time when the CF velocity at the cylin- 
drical wall exceeds a constant critical value. Fig. [j] il- 
lustrates a typical final configuration for such a cluster 
which includes only one species of vortex lines. They 
are introduced at a fixed nucleation center, once the CF 
velocity reaches the critical value. 

To simplify matters we may here think of the nucle- 
ation center in the form of a line source parallel to the 
rotation axis, located close to the cylindrical wall of the 
container From the nucleation center the newly cre- 
ated rectilinear vortex line moves to the edge of the vor- 
tex cluster under the action of the background potential 
of uniform rotation along a spiral trajectory, which de- 
pends on the values of the mutual friction parameters b 
and b' . The line is incorporated into the cluster by push- 
ing via the repulsive interactions, which leads to plastic 
deformations in the existing array. The final structure 
of the array resembles in the central region a hexagonal 
lattice with a large number of dislocations. The lines in 
the periphery form circular rings which are separated by 
a wide transition region from the central crystal. 

A second variation of the same calculation, but now for 
a two-component vortex array, is shown in Fig. g. Again 
the final state is shown, after DQV lines have been ac- 
cumulated one by one into an existing cluster of SQV 
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FIG. 1. Computer simulation of vortex array formation. 
Here Ni = 127 rectilinear SQV lines have been accumulated 
one by one to form a cluster, by increasing Q slowly in the 
superfluid state. The equation of motion ([ll]) with 6' = 
has been used. The nucleation center, from where the lines 
are injected, is located at the cylinder wall at x = R (not 
shown) and y — 0. The final configuration of the array in the 
cylindrical container is shown. The dashed circle marks the 
edge of the cluster at r = R v . In the center the lines form 
a 2D hexagonal crystal with a large number of dislocations, 
whereas towards the perimeter the configuration deforms to 
circular concentric rings. 



lines. Despite large plastic deformations in the cluster, 
the injected DQV vortex lines form a compact drop. We 
observe neither a periodic (chemical) or irregular struc- 
ture of single DQV lines inside the SQV array. Instead 
this simulation procedure gives a glassy vortex agglom- 
erate with hexagonal coordination, which is consistent 
with the experimental situation in 3 He-A (as will be seen 
below by comparing eg. with Fig. 14). 

To explain these results let us first consider the energy 
of mixing in different periodic vortex structures, with two 
types of vortices in the unit cell. Fig. |l| shows the differ- 
ence between the free energy (per vortex line) of some pe- 
riodic structures (with up to 9 vortices per unit cell) and 
the sum of the energies of two separate perfect hexagonal 
crystals formed from SQV and DQV lines. The values in 
Fig. ^ have been obtained using the planar summation 
method from Ref. |l4| ]. We see that for all these struc- 
tures the energy of mixing is positive - the lowest energy 
corresponds to two simple hexagonal separated lattices. 
This fact explains the tendency to form a compact pocket 
of DQV lines inserted into the SQV array. 

The second feature of the simulation is that the equa- 
tion of motion (Eq. (|ll|)) does not include any thermal 
annealing process. This, in fact, corresponds to the situ- 
ation in superfluid 3 He: Thermal fluctuations are weak, 
because of the very low temperatures compared to high 
energy barriers, associated with variations in the super- 
fluid velocity field. The above arrays are therefore not 




FIG. 2. Simulation calculation of an array with 100 SQV 
and 50 DQV lines. The DQV lines have been accumulated 
one by one to an existing equilibrium array of SQV lines, sim- 
ilar to the experiment in Fig. [u| The nucleation center, from 
where the DQV lines are injected, is located on the cylindri- 
cal wall (x = R, y = 0) and is marked with a cross in the 
inset where the final configuration of the array at 0.49 rad/s 
is shown. With this number of vortices the array would cor- 
respond to equilibrium circulation at 0.34 rad/s. The thin 
curves in the main panel represent the radial locations of the 
SQV and DQV lines in the array: N a (r) = Ylji x a,j + Vaj)- 
They illustrate how the locations of the vortex lines concen- 
trate with increasing radius more and more on concentric 
rings, although the nearest-neighbor coordination is still tri- 
angular, as seen in the inset. The thick curves correspond to 
deceleration records (fi : 0.49 rad/s — > 0), similar to those 
measured in the experiment (Fig. Nj), and give the average 
radial distributions Ni(r 2 ) and Afofr ). 



expected to have the optimal structure of the global free 
energy minimum. 



2. Equilibrium arrays 

Experimentally the above simulations correspond to a 
situation where the container with superfluid 3 He-A is 
slowly accelerated to rotation at constant temperature 
below T c . Cooling the container through T c in rotation 
at constant Q is expected to produce a different situation, 
the global free-energy minimum. To find this state, we 
now compare the free energies of different vortex array 
configurations. First we note that for evaluating integrals 
such as those in Eq. ( |lO| ) it is useful to introduce a vortex- 
core cutoff at short distances such that close to the core 
of the jth vortex we replace 
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FIG. 3. Energy of mixing SE of different periodic 2D 
vortex line lattices, formed from SQV and DQV lines. The 
zero level corresponds to the energy of two separated hexag- 
onal SQV and DQV lattices. The unit cell includes Ni SQV 
and N2 DQV lines. The horizontal axis gives the fraction 
p = Ni/Nc, where N c = Ni + N2. The energy of mixing 
5E is expressed in units of p s Kq/(2-kN c ). The lattices with 
N c = 2 and N c — 4 have a square unit cell, while all other 
lattices are obtained by inserting the additional vortex lines 
into a hexagonal cell. 



where | r — Vj | < ^ . Using Eq. (|io| ) and performing 
integration by part, while neglecting all terms of order 
(Ca/^a) 2 since £ a <§; r a is assumed, we can rewrite Eq. ([]]) 
in the form 
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Next it becomes useful to introduce a vortex lattice 
energy which provides a more sensitive measure of the 
configuration dependent energy differences. Let us con- 
struct the same free energy as in Eq. (Q) but with the 
superfluid velocity written in the continuum approxima- 
tion with v s = v s o, where v s o is defined in Eq. (^|). The 
free energy in the continuum approximation is then (cf. 
Eq. (I)) ' 



£{v s o} = ^ / (vf - 2v s0 v rot )d 2 r 
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The difference of the two energies is the Madelung energy, 



E M = E{v s } - E{v s0 } 



(15) 



The Madelung energy represents the vortex lattice con- 
tribution and includes the interactions with the image 
vortices. If one neglects the difference v s — v s0 in the 
vortex-free region, 
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Within the logarithmic approximation it reduces to the 
last term in Eq. (||). 

An important feature is the interaction with walls. By 
moving the wall further from the vortex cluster its influ- 
ence is reduced and the vortex cluster starts to approach 
"ideal matter" with more universal properties, similar to 
point charges in an ion trap. Consider the outermost ring 
of vortex lines in an ideal coaxial cluster. Nonuniformi- 
ties in the velocity field around the cluster, which arise 
from the discrete nature of the vortex lines in the outer 
ring, decay as ~ exp(— 2 / n5r/r a ) with distance Sr from 
the outer boundary ]l5| . Conversely, the interaction of 
the outermost vortex ring with its reflection is of order 
exp(— 4:Ttd/r a ), where d is the distance between the ring 
and the container wall. To reduce this effect, we take the 
container radius R to be clearly larger than R v in our 
simulation calculations. We then obtain an estimate for 
the Madelung energy which does not depend on the inter- 
action of the particular vortex array with the container 
wall. 

Let us consider an array consisting of only one type of 
vortex lines (e.g. a). We expect that ideally the mini- 
mum energy at large N a is reached with a configuration 
with hexagonal structure in the center of the cluster, cir- 
cular rings at the outer boundary, and an intermediate ir- 
regular region between these two regimes. In Eq. ( 16) v s o 
is the average superfluid velocity within the cluster, equal 
to the velocity of rotation (Eq. (0)), and thus v s — v iS o 
represents the deviation from the average in the rotating 
coordinate system. In an ideal periodic structure this ve- 
locity is also a periodic function of the coordinates and 
the integration in Eq. ([l6]) is reduced to one over a single 
unit cell: 
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As a first approximation we may replace the hexagonal 
lattice cell with a cylindrical Wigner-Seitz cell of the same 
area. The area per vortex is b 2 \/3/2 = nr 2 , which gives 
b a = 1.9046r a as the relation between the lattice con- 
stants b a of a 2D hexagonal lattice and r a of a cylindrical 
WS lattice. We refer to this approximation as the WS 
cylinder. The periodicity of the lattice requires that the 
component of v s — v s q, which is normal to the border 
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of the unit cell, vanishes at the border. In the rotat- 
ing coordinate frame the superfluid velocity has only an 
azimuthal component and we may write 



v s0 
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The Madelung energy can then be written as the energy 
of a single vortex 
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where olm will be referred to as the Madelung con- 
stant. For the WS-cylinder approximation in Eq. (18), 
a M = 1/2. Comparing Eqs. @ and ( |l9| ) we obtain 
the expression for the energy density e a through the 
Madelung energy as in Eq. (Q). The Madelung constant 
«m and the Madelung energy can be calculated exactly 
using the method of planar summation as in Ref. [fl4|| . 
For an ideal hexagonal structure it yields aM = 0.49877. 
This value is very close to that of the WS-cylinder ap- 
proximation: aM = 1/2. 

The Madelung energy in Eq. ( |l5| ) for large, but fi- 
nite N a can be expanded in powers of the small variable 
b a /R v : 
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Here the first term within the square brackets is the vol- 
ume contribution, and the term proportional to (3 is the 
surface energy per unit vortex length 
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The number of vortex lines in the outermost ring is 
27r_R„/6 Q , where the intervortex distance b a is roughly 
equal to the lattice constant of the ideal hexagonal struc- 
ture: b a w 1.90r a . One can check that the factor (3 is 
independent of the type of vortex as long as the inter- 
action of the vortex lines with the cylinder wall can be 
neglected. 

To determine the surface energy factor /?, we calculate 
Em numerically for different values of N a . This requires 
minimization of the free energy in Eq. ( |l3| ) for differ- 
ent values of N a . Here we use the simulated annealing 
method from Ref. where it was applied to calculate 
the energy of a large 3D ion cluster in an ion trap. First 
the vortex lines with a fictitious mass are heated to some 
"temperature" which is modeled by a random force. The 
cooling is modeled by inserting a small viscous term in 
a "Newtonian equation of motion" . The heating-cooling 
cycles are repeated, to allow the array to converge, as a 
function of time, towards some minimum-energy config- 
uration with a fixed number of vortex lines. Eq. ( |l3| ) 
is then used to calculate the free energy of the final 



converged optimum configuration. Its Madelung energy 
is obtained from Eq. (|l5|). Note that plastic deforma- 
tions play an important role in making relaxation pos- 
sible towards the energy minimum. The (logarithmic) 
2D Coulomb interaction can be regarded as a soft-core 
potential, in contrast to a hard-sphere interaction of eg. 
the Lennard- Jones type. This simplifies the relaxation of 
the vortex array towards an optimal structure. Moreover, 
in contrast to 3D Coulomb systems, where many crystal 
structures (bcc, fee, hep) have nearly the same energy, 
the energy of a 2D hexagonal crystal is lower than the 
energy of most other structures. 

Fig. H shows the Madelung energy Em /N\ per vortex 
as a function of n/ R v = 1/y/Nx, for arrays consisting of 
SQV lines. We start the simulation from an ideal hexago- 
nal lattice within a circle r < R v and after many (~ 10 2 ) 
heating-cooling cycles we arrive at the structure which 
corresponds to a stable minimum of Eq. (|l^). By re- 
peating this procedure for different values of R v (or N\ 
we obtain j3 in Eq. ( po|) from the fitted slope in Fig. 
(3 = 0.010. The calculated optimal structure of these vor- 
tex clusters (Fig. ||) approaches that of an ideal hexagonal 
lattice in the center, but then distorts to a quasi-circular 
structure towards the outer boundary, with a transition 
region between these two extremes. 

Similar simulations for the equilibrium state of a clus- 
ter with two coexisting vortex species are more delicate. 
First, we have already noted that the energy of a mixture 
of SQV and DQV lines is higher than that of two segre- 
gated domains. Secondly, it follows from Eq. ( ^l| ) that 
the surface energy of a DQV cluster is 2y/2 times larger 
than that of a SQV cluster with the same value of R v , i.e., 
02/ci = 2-\/2 because k 2 = 2«i and r 2 — v2ri- Hence, 
the only configurations, which were found to converge 
towards a stable minimum-energy state, are coaxial clus- 
ters with the DQV lines in the center, surrounded by an 
outer annulus of SQV lines. In Fig. g the Madelung en- 
ergy Em /N per circulation quantum of such clusters has 
been plotted vs. n/R v = 1/VN, where N = Ni + 2N 2 
is the total number of circulation quanta. We take here 
Ni = 2N2 and thus the interface separating the two types 
of vortex lines lies at r « R v /V2. The optimal structure 
(Fig. 0) turns out to consist of a nearly ideal hexagonal 
lattice of DQV lines, which extends all the way from the 
center to the inner interface. At the inner interface the 
mismatch in periodicity between the two vortex species, 
b 2 = y/2h, creates an abrupt break, but the smaller SQV 
lines very effectively compensate for the faults on the 
boundary of the central hexagonal DQV crystal. Also 
the SQV lines fulfill the boundary condition and on ap- 
proaching the outer boundary of the cluster a perfect 
circular ring structure is established. 

Simulation attempts on the inverse initial state, 
namely on segregated coaxial domains but with the SQV 
lines in the center, were not successful. Such calculations 
are slow in convergence since they start to drift towards 
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FIG. 4. Madelung energy Em for an equilibrium cluster 
of SQV lines, plotted as a function of ri/R v — 1 /yNi (where 
Ni — 55, 85, 151, 199, and 253). The calculations represent 
the minimum value of Em/Ni (in units of p s kq/2tt) after a 
large number of simulated annealing cycles. The container 
radius is set at R = 40ri while the core radius is fixed to 
£i/n = 0.218. The latter corresponds to SQV lines with 
£i = 32/im at = 0.5 rad/s. The dashed line represents the 
linear dependence of Eq. (EOh and the fitted slope gives for j3 
the value : (3 = 0.010. The extrapolation Em/Nx = 0.512 to 
an infinite vortex array equals that which one obtains from 
Eq. (Jl9|) for an ideal hexagonal lattice. 



the ideal structure where the SQV lines are in the outer 
periphery. The same applies to asymmetric starting con- 
figurations, such as the one where the cluster is divided 
along its diameter into equal halves with only SQV and 
DQV lines in each. This arrangement does not minimize 
the outer surface energy or allow for the presence of a 
symmetric hexagonal crystal in the center of the array. 

We can use the results in Fig. || to estimate the sur- 
face energy of the inner interface a±2- The dashed line 
shows a fit to the calculations: E M (DQV + SQV)/N = 
0.7294 + 0.0462 n/R v (in units of p s n 2 /2ir). The same 
fit to the SQV clusters in Fig. | gives E M (SQV)/N = 
0.5122 + 0.0348 n/R v . Since the outer shell in both 
cases consists of SQV lines only, we can associate the 
difference in slopes with the surface energy of the in- 
ner interface between the SQV and DQV domains. This 
gives 6E12/N = 0.012 ri/R v . Taking into account the 
length 27ri?2 ~ 2ttR v /^/2 of the interface, we obtain 
f 12/f 1 ~ 0.48. Remembering that a^/ci = 2.83, we find 
that o"2 > <7i + cri t 2- Thus the contact angle between the 
SQV domain and the vortex-free CF annulus has to be 
zero: the SQV domain "wets" completely the vortex-free 
region. The most preferable configuration is therefore 
one where the DQV domain is completely surrounded by 



FIG. 5. Calculated minimum-energy structure of a SQV 
cluster with Ni = 253 lines. This is one of the examples in 
Fig. ^. The dashed circle marks the edge of the cluster at 
r — R v . The container wall is further out at R = 2.5R V . 
This is thus a metastable state with less than the equilibrium 
number of vortex lines, but with the cluster itself at equilib- 
rium density. All calculated clusters in Fig. ^ are of this type, 
to reduce the influence of the cylinder wall on the simulation 
calculations. 



SQV lines. 

Interestingly the reduced surface energy with SQV 
lines makes the coexistence of DQV and SQV lines pos- 
sible even in a region of the SQV - DQV phase diagram 
where the global energy minimum corresponds to that 
of DQV lines only. The energy of a cluster consisting of 
DQV lines only is given by Eq. (f20|), 

(2 ) PsKjRl ( r 2 \ 
E M = ~A — 2~ m 7 a M + 2nR v <7 2 , 

while the energy of a cluster with a central DQV domain 
surrounded by SQV lines is 

(2i) Psi^Rl ( r 2 \ 

H — * m 7 a M 

+2irR v <Ji + 27ri?2Ci2 ■ 

(21) 

E M as a function of R 2 has a minimum at R 2 = for 
il < fT, when a SQV cluster is more favorable. However, 
due to a finite a± 2 , there is still a minimum at R 2 = even 
when > Q*, i.e. when a DQV cluster is more favorable. 
In this second case another minimum appears at R 2 = R v 
which corresponds to a DQV cluster surrounded by a thin 
surface layer of SQV lines. Indeed, the energy E y M ' is 

a decreasing function of R 2 in the region R 2 ~ R v . For 

(2) 

R 2 = R v it is smaller than E y M ' because a 2 > o\ + o\ 2 . 
The difference SE = E^ 1] - E { B for ft > Q* is shown in 
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FIG. 6. Calculated Madelung energy per circulation quan- 
tum, Em /N, for segregated coaxial minimum-energy clusters, 
composed of SQV and DQV lines, plotted as a function of 
the inverse of the normalized cluster radius: r\/R v = l/y/~N. 
In the minimum energy configuration the DQV lines are in 
the center. Here the total circulation of the DQV lines 
has been taken to be equal to that of the SQV lines (i.e. 
iVi = 2N2 — N/2). Em is expressed in units of p s kq/2tt and 
is calculated for five cluster radii N2 = 55, 85, 121, 151, and 
199. The parameters are £i/n = 0.218, £ a /ri = 0.236, which 
correspond to £1 = 32 pm and £2 = 49 jim. at fi = 0.5 rad/s. 
The container radius is set to R — 40 n = 5.8 mm. The arrow 
indicates the sum of the Madelung energies in Eq. ( ]2o| ) for in- 
finite SQV and DQV lattices, (E M i + E M 2)/N in the limit 
N — > 00, when the surface energy can be neglected. The solid 
fine has the same slope as the line in Fig. Q and represents the 
surface energy of the outer boundary of a cluster with SQV 
lines. The dashed line is a linear fit through the data points. 



Fig. H as a function of R2/ R v - The two local minima at 
i? 2 = and i? 2 ~ Rv have the same energy when 



fit = fi* 



1 



8irf3rl a 12 



Rvh o\ 



(22) 



Here fi* is determined by Eq. (g) where, however, e 
should be replaced with the more accurate expression 
e 2a M f an idgaj vortex lattice, fit is only slightly higher 
than fi* because of the large value of R v jr\. 

Eq. ( P4 ) states the condition for the global equilibrium 
between SQV and DQV lines. However, because of the 
minimum at R « R v , the model allows for the coexistence 
of the two vortex species: A "monolayer" in the form 
of an outer circle of SQV lines could surround a DQV 
cluster, when fi > fi*. Evidently, our idealized model 
only applies if the vortex layer is at least a few inter- 
vortex distances wide. With this warning we note that 
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FIG. 7. Calculated minimum-energy structure of a cluster 
with N2 = 199 and Ni — 2N2 lines. This is one of the exam- 
ples from Fig. [| The dashed circle denotes the edge of the 
cluster at r = R v while the container wall is at R — 1AR V . 



the qualitative behavior of SE is shown in Fig. || by the 
solid line, when an outer boundary layer of SQV lines is 
present. 

Summarizing our analysis we conclude that the mini- 
mum energy configuration for coexisting SQV and DQV 
lines in a rotating cylinder can only be established by 
considering discrete vortex arrays. The optimal configu- 
ration (1) consists of segregated domains, (2) has cylin- 
drical symmetry (3) with the SQV lines forming the outer 
boundary of the cluster. 



IV. EXPERIMENTAL METHOD 

The structure of the vortex cluster in the rotating con- 
tainer is probed with cw NMR. It allows one to determine 
the average radial composition, ie. Ni(r 2 ) and A^r 2 ). 
The NMR method is based on a measurement of the ab- 
sorption spectrum to which SQV and DQV lines each 
contribute a characteristic satellite peak. These satel- 
lites have frequency shifts which are different from each 
other and from the bulk liquid NMR line |ll| . Both the 
peak height and the integrated intensity of the satellites 
are proportional to the total number iV a of vortex lines of 
the species (a) and can be calibrated to give the absolute 
value of N a . 

To measure the composition of a vortex cluster we 
monitor the peak height of a satellite as a function of the 
rotation velocity fi during slow deceleration, when the 
cluster expands and the outermost vortex lines annihi- 
late one by one at the cylinder wall [|l7j . The peak height 
as a function of fi gives the number of the annihilating 
vortex lines as a function of r 2 . We can thus reconstruct 
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FIG. 8. Free energy difference 5E between a DQV cluster 
with radius R2 and one which has an outer surface layer of 
SQV lines in the region R2 < r < R v , when SI > SI*. Within 
the idealized model of Eq. , there exists a step-like energy 
reduction by 2nR v {a2 — o~\ — (T12) (dotted line), when a thin 
layer of SQV lines is placed on the outer surface of a DQV 
cluster. Decreasing the radius R2 of the DQV domain further 
results in an increase in the free energy ( dashed curve ) until it 
reaches a maximum followed by another minimum at R2 = 0. 
The solid curve indicates schematically the dependence with 
a SQV surface layer of finite width. 



the cumulative number of both types of vortices N a {r 2 ), 
which gives the number of vortex lines of type (a) inside 
a radius r. Indeed, for reconstructing the radial distri- 
butions only one of the peak height dependences N a (r 2 ) 
needs to be measured. The requirement of solid-body 
rotation of the superfluid component within the vortex 
cluster gives the relation: Ni(r) + 2N2(r) — 2irr 2 Vl/ k$. 

The measurements have been performed with two 
cylindrical containers, one fabricated from epoxy and the 
other from fused quartz Both arc right circular 

cylinders 7 mm long and R = 2.5 mm in radius. They 
are closed off, except for an orifice of 0.5 mm in diame- 
ter, located in the center of the flat bottom plate. The 
orifice connects to a long tubular channel which provides 
via its liquid 3 He column the thermal contact to the re- 
frigerator. The 3 He-A temperature in the container is 
determined from the calibrated temperature dependent 
frequency shift of the bulk liquid NMR absorption peak. 

Two controlled ways exist for preparing a vortex clus- 
ter: 1) The sequential addition of vortex lines one by one 
to the cluster by accelerating £1 slowly (in the presence 
of a constant critical velocity), and 2) by cooling in ro- 
tation through T c . One might expect in the latter case 
that a sufficiently slow cool down through T c , one which 
approaches the adiabatic limit, will lead to an equilib- 
rium state, where the number and type of vortex lines as 
well as the configuration of the cluster adjust such that 
the total energy becomes minimized. This has been ex- 
perimentally checked with respect to the type of lines || 
and the number of lines |L7| . 

In 3 He-A the transition between the SQV and DQV 
lines occurs in the presence of an applied magnetic field 



pT| as a function of £1 g]. Experimentally we find that 
there exists a region of coexistence around this transi- 
tion, ie. the transition has a width as a function of ft 
within which both species of vortex lines are present si- 
multaneously, but where their relative amounts depend 
on £1. Ideally the transition should be sharp and there- 
fore our experimental situation cannot represent the ho- 
mogeneous adiabatic limit. Nevertheless, in the coexis- 
tence regime the energy difference between the two vortex 
species is small, and, in principle, it can be smaller than 
the energies associated with the vortex cluster configu- 
ration. Therefore, in addition to trying to identify what 
type of inhomogeneity controls the width of the transi- 
tion, we may investigate the vortex-array configurations 
in the coexistence regime and ask whether these energy 
contributions are minimized in the transition. 

The sequentially formed coexistence clusters are more 
obvious. Only one type of cluster is possible: Due to 
their much higher critical velocity SQV lines can only be 
formed by cooling through T c in rotation. Thus by cool- 
ing slowly through T c at low £1 an equilibrium cluster of 
SQV lines is formed, to which DQV lines can be added 
sequentially by increasing £1 slowly at constant tempera- 
ture below T c . 

The first observation from our measurements was that 
the radial distributions of the two vortex species turned 
out to be different in the equilibrium clusters and the se- 
quentially constructed states. This means that in super- 
fluid 3 He a common cluster configuration is not realized, 
but metastable states are long lived and the structure of 
a vortex cluster depends on how it has been prepared. 



V. MEASUREMENTS 
A. Coexistence regimes of SQV and DQV 

1. Dependence on rotation while cooling through T c 

When the rotating sample is cooled in the NMR field 
through T c , an equilibrium vortex state with only SQV 
lines is formed at low rotation, £1 < 0.6 0. Above this 
velocity also DQV lines are created. Their relative num- 
ber increases with £1 until above 3 rad/s no more SQV 
lines are observed. The number of DQV lines is shown 
in the coexistence regime as a function of £1 in Fig. [{j, 
where each data point (o) represents I v /Itot, the inte- 
grated NMR absorption I v of the DQV satellite normal- 
ized to the total absorption I tot in the NMR spectrum. 

The normalized absorption I v (fl)/I to t can be trans- 
lated to vortex line numbers A^(n) by comparing to an 
equilibrium state with only DQV lines. The equilibrium 
number of vortex lines (of type a only) is given by 



A a(eg) (£l) = tt(R - d a(eq) ) 2 — cx -y 



2 ™ cx l 4± , (23) 

Kn 
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FIG. 9. Coexistence of SQV and DQV lines in the 
high-field (dipole-unlocked regime as a function of ro- 

tation velocity S2. (o): Normalized DQV satellite inten- 
sity I v /Itot plotted as a function of the rotation velocity at 
which the sample was cooled through T c . On the right ver- 
tical axis the intensity has been expressed in terms of line 
numbers, N 2 (Q), by means of the calibration, N 2 ( eq ). At 
low Q only SQV lines are formed and N2(fl < 0.6 rad/s) 
= 0. Above 0.6 rad/s the coexistence regime starts, where 
the fraction of the total circulation in DQV lines increases 
with Q, until well above 3 rad/s no more SQV lines are 
formed. As a guide for the eye, the solid curve represents 
h/Itot = 0.033(S1 - 0.6) + 0.0084(S1 - 0.6) 3 . The cool-down 
rate at T c has been kept at \dT/dt\T c < 5 /iK/min, ex- 
cept for the data point marked with an open square (□), 
for which the cooling rate was (dT/dt)T c = —18 ^K/min. 
(•): Calibration measurement of I v /Itot with the equilib- 
rium number of DQV lines, Nweq) ■ The fitted result gives 
I v (Sl)/Itot = 2.74 • 10" 4 N 2 ( eq )(Q) « 0.0806 fi (1 - O.l/v 7 ^) 
(fl in rad/s). The inset shows the radius of the central DQV 
cluster in relative units in the coexistence regime, together 
with a guide for the eye: R 2 /{R-d 1{eq) ) = 0.618 (fi-0.6) ' 41 . 
Other conditions: P = 29.3 bar, H = 9.91 mT, the satellite 
intensities have been measured at T = 0.81 T c . 



where the width d rwj of the equilibrium CF annulus 
is defined in Eq. (|9|). The calibration measurement, 
iV 2 ( eg ) (f2), is plotted (•) in Fig. ^|. In practice, the equilib- 
rium number of vortex lines is equal to that at the annihi- 
lation threshold [[17j . In Fig. [)] the data for N 2 ( eq \ (fi) were 
measured by first slowly increasing f2 from zero to some 
high velocity f2 > 3.5 rad/s at T = 0.81 T c and by then 
decelerating sequentially to lower SI values and recording 
there at constant SI the respective value of I V (U) / Itot- 
The difference between the two curves in Fig. || gives 
the number of SQV lines in the metastable coexistence 
regime: N^Q) = 2 [N 2{eq) (Q) - JV 2 (n)]. 

The critical rotation velocity S7* for the transition from 
SQV to DQV lines lies thus according to Fig. |^ within 
the interval 0.6 - 3 rad/s. An estimate for SI* can be 
worked out from Eq. (g) using the normalized satellite 
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FIG. 10. Coexistence of SQV and DQV lines as a function 
of magnetic field H below 0.6 rad/s: Normalized DQV satel- 
lite intensity I v / hot plotted vs the field H, which was applied 
during cooling through T c at 0.5 rad/s. On the right verti- 
cal axis the NMR absorption intensity has been converted to 
DQV vortex numbers, N 2 (H), using the calibration in Fig. []. 
The horizontal dashed line represents the equilibrium number 
N 2 ( eq ) of DQV lines. The shaded region marks the width of 
the coexistence regime around the critical field, H c m 0.6 mT. 
In this measurement the field was increased to the NMR value 
9.91 mT in the superfluid state below T c , while the sample 
was cooling to 0.81 T c , where the intensity I v /Itot was mea- 
sured. Other conditions were the same as in Fig. In the 
transition through T c the equilibrium value of total circula- 
tion is created and thus at high field H > H c , where N 2 — > 
in accordance with the results in Fig. ^, the deficit from N2( eq ) 
represents the SQV lines. 



intensity I v /Itot- It measures the effective vortex-core 
area compared to that of the WS cell [O] , or 



A. 

Hot 



(24) 



For the sake of comparison we neglect the temperature 
dependence of the core radii (jll| and use the best val- 
ues for the satellite intensities, which come from mea- 
surements at 0.50 T c (11.7 mT, and 33.9 bar): In the 
equilibrium state the coefficient in front of ft in Eq. (|24| ) 
was found to be 0.094 for SQV and 0.114 for DQV lines. 
These values give for the soft core radii £i w 32 /im and 
£2 ~ 49 fim. From Eq. (||) we now get CI* w 2.8 rad/s, in 
qualitative agreement with the experiment. The estimate 
for f2* in Eq. (^) is based on the differences in quantiza- 
tion number, core size, and the crude approximation for 
vortex energy in Eq. (|t]) . A more sophisticated numerical 
calculation in Ref. || places ST in the interval 0.54 - 1.4 
rad/s. 
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2. Dependence on magnetic field while cooling through T c 

In addition to fl, the second important variable ol the 
vortex phase diagram is the applied magnetic field H 
fL9f : When it exceeds the equivalent of the spin-orbit 
coupling, as is the case in our NMR measurements, the 
SQV and DQV structures display well-separated vortex 
cores, with r a £ a jy]]. At low fields the spin-orbit in- 
teraction dominates and these vortex structures are not 
stable. Instead here vorticity appears in the form of ex- 
tended periodic vortex textures with no core or a much 
expanded core, such that r a ~ £ . The dominant low- 
field vortex texture has a square lattice and a unit cell 
with four quanta of circulation J2|,^| . These vortex tex- 
tures are separated from the high-field structures by a 
transition of first order. 

The transition between the low and high field regimes 
as a function of H is illustrated by the measurements 
in Fig. [To]. Our NMR is performed at a fixed frequency 
in the high field limit. Therefore the measurements in 
Fig. [To] have to follow a particular routine. First, the 
container is cooled through T c at 0.5 rad/s, with the 
field H adjusted to the value given on the horizontal 
axis. Next, below T c the field is increased to the NMR 
value of 9.91 mT. During the field sweep the low field 
vortex texture undergoes a topological transition to the 
high-field (dipole-unlocked) DQV structure, with well- 
separated vortex cores. This transition is not revealed in 
our measurements, if the field sweep is conducted at high 
temperatures (T > 0.8 T c ). However, the SQV has differ- 
ent topology [ ]ll| and any vortices of this kind will remain 
singly quantized during the field sweep. After the field 
sweep, the container is cooled to lower temperatures, to 
improve the resolution in the NMR measurement of the 
vortex satellites, and the DQV satellite intensity is mea- 
sured. 

As shown in Fig. [lO|, we find a critical field H c « 0.6 
mT. Above H c the SQV lines dominate, as required when 
fl < 0.6 rad/s. Below H c only DQV lines are found, 
which represent the transformed low-field vortex textures 
after the field sweep. Centered around H c there is a 
coexistence regime with a width of 0.5 mT, where both 
types of vortex lines are present. Within the coexistence 
regime, the relative abundance of SQV lines increases 
continuously with increasing H , similar to the situation 
in Fig. |asa function of f2. 

The same behavior is verified in Fig. |ll| at a higher 
rotation velocity of 1 rad/s. On comparing Figs. [lO] and 
[TT| it is noted that both H c and the width of the coexis- 
tence regime increase with Q |2| . In the measurements of 
Fig. [ll] both the DQV and SQV satellite intensities were 
determined separately. On the low field side only DQV 
lines are found. On the high field side both types of lines 
are formed, with a 4 times larger share of the circulation 
carried by the SQV lines. This is the same ratio as was 
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FIG. 11. Coexistence of SQV and DQV lines as a function 
of magnetic field above 0.6 rad/s: The normalized intensities 
Iv/hot of the DQV satellite (o; left vertical axis) and I s /Itot of 
the SQV satellite (•; right vertical axis) have been plotted as 
a function of the magnetic field H , which was applied during 
cooling through T c at 1 rad/s. The horizontal dashed line on 
the left represents the equilibrium value expected for DQV 
lines, after a superfluid transition at low fields (H < H c ) and 
a subsequent increase of the field to 9.91 mT. On the right 
(H > H c ), the dashed lines give the mean values for the SQV 
and DQV lines: Ni « 440 and 7V 2 w 48. The shaded region 
marks the width of the transition regime around the critical 
field, H c w 1.2 mT. The coexistence regime includes both 
the shaded and the high field regions, H > H c . The same 
conditions apply as in Fig. ^| including the calibration for the 
DQV satellite intensity. The calibration for the SQV satellite 
intensity is N 1(eq) {Q) tx I 3 (Q.)/I tot « 0.43 0,(1- 0.089/^0) 
(fi in rad/s). 



already seen in Fig. ^ at 1 rad/s. 

In a magnetic field the superfluid transition splits into 
T c i and T C 2- The former is the transition from the normal 
to the Ai phase and the latter from Ai to the A2 phase, 
which is the regular 3 He-A. The thermal width of the Ai 
region is narrow, approximately 0.6 /iK at 10 mT. We 
might wonder what effect the intermediate Ai phase has 
on the DQV to SQV ratio above H c in Figs. g - [0]: Are 
more changes to be expected at fields much larger than 
those shown in Figs. [l(] and [Tl"| ? We find that at least 
a cool down in 60 mT field at 1.4 rad/s gives the same 
result as at 10 mT. 



3. Dependence on cooling rate at T c 

Surprisingly in Fig. |^ one data point deviates from the 
typical behavior (denoted with (□) at 1.6 rad/s). The 
only distinction that can be appended to its measure- 
ment, compared to other data points at the same f2, is a 
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FIG. 12. Influence of cooling rate at T c on the relative 
number of SQV and DQV lines. The left vertical axis gives 
the normalized DQV satellite intensity I v /Itot as a func- 
tion of cooling rate during the transition through T c at 1.4 
rad/s. The satellite intensities are measured at (0.49 - 
0.51) T c , where the calibration for the DQV line satellite is 
N 2(eq) (n) oc I v (n)/I tot « 0.078 fi (l-0.ll/-/n) (n in rad/s). 
In this expression I v includes only the primary satellite peak 
(the secondary peak, which is also visible at these low tem- 
peratures, is not included in the absorption). On the right 
vertical axis the calibration has been used to list the fraction 
of the total circulation in DQV lines, N 2 /N 2 ^ q - j , while the 
rest, 1 — N 2 /N 2 ( eq - ) , is carried by SQV lines. 



larger cooling rate during the superfluid transition: The 
five data points measured at cooling rates \dT/dt\r a = 1 
- 5 /iK/min give for the relative amount of DQV lines 
N^/N^Uq) ~ 0.38 while the anomalous point at a cooling 
rate of 18.5 /dC/min falls on the equilibrium curve N 2 ( eq ), 
with no SQV lines detected at all. 

Additional evidence is shown in Fig. [l2]on how a faster 
cooling rate at T c shifts vorticity from SQV to DQV lines. 
In these measurements a fused quartz container was used 
in which 3 He-A could be supercooled to below 0.49 T c . 
This feature improves the measuring resolution since at 
lower temperatures the different vortex satellites are sep- 
arated by larger shifts from each other in the absorption 
spectrum |[(|. When the cooling rate in Fig. [l2] at 1.4 
rad/s reaches 5 /iK/min a change over takes place from 
SQV to DQV lines: at lower rates only one fourth of 
the total circulation is in DQV lines (which is the same 
fraction as in Fig. |^), while above 6 /iK/min, only few 
SQV lines are formed. Similar measurements at 0.6 rad/s 
with cooling rates up to 16 /iK/min did not yet produce 
a DQV satellite of measurable height. 

The quartz container in the measurements of Fig. [l^ 
was thermally connected to the refrigerator via an ori- 
fice of 0.5 mm in diameter. Later measurements with 
another quartz cylinder with a small and long channel of 




FIG. 13. Radial distributions N^r 2 ) and N 2 (r 2 ) after 
cooling through T c at 1 rad/s and 1.1 mT. The deceleration 
from 1 rad/s to zero is performed at 0.81 T c and 9.91 mT. 
Other conditions are the same as in Fig. jw| N 2 ( eq ) shows for 
comparison the deceleration for an array consisting of only 
DQV lines. It is prepared by cooling through T c at > 1.5 
rad/s in zero field and by then sweeping the field up to 9.91 
mT at 0.81 T c . The final part of the deceleration below 0.7 
rad/s is shown here. As demonstrated in Ref. J17J, during 
a slow deceleration (d£l/dt = —1.5 ■ 10~ 3 rad/s - ^the total 
circulation at any value of il corresponds roughly to that in 
the equilibrium state. The inset illustrates schematically a 
transverse cross section through the rotating container with 
the coaxial array configuration: the SQV lines (small dots) 
are in the outer annulus and the DQV lines (large dots) in 
the center. 



0.4 mm diameter displayed appreciably smaller fractions 
of SQV lines, also in the regime f2 < 0.6 rad/s where 
N 2 = in Fig. |. The narrow and long channel increases 
the thermal resistance between the sample container and 
the refrigerator, and leads to rapid irreversible normal-to- 
superfluid transitions, where the transition velocity can- 
not reliably be controlled externally. This suggests that a 
large cooling rate and nonequilibrium conditions enhance 
the formation of DQV lines at the expense of SQV lines. 

There is one data point also in Fig. at 6 mT with a 
larger fraction of DQV lines than what one would expect 
on the basis of the general behavior. In the measurements 
as a function field (Figs. [Io| and [yj) the cooling rate at 
T c cannot be monitored accurately, since the NMR mea- 
surement is inoperative at that stage. Therefore a cool- 
ing routine was worked out which should have yielded 
rates of less than 5 /iK/min. However, without the NMR 
record it is not possible to state whether a mistake oc- 
curred in the cooling routine of the anomalous data point 
in Fig. fiR 
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B. Vortex cluster in coexistence regime 

The radial composition of a cluster with both SQV 
and DQV lines can be determined by recording the am- 
plitudes of the satellites as a function of £1 while the rota- 
tion is slowly decelerated. In Fig. [l3] the result is shown 
for an array which has been created by cooling through 
T c at 1 rad/s. Here the distributions of the two vortex 
species are segregated: The SQV lines annihilate first 
and the DQV lines later. This behavior corresponds to 
the coaxial array structure, shown in the inset of Fig. |l3| , 
which consists of a central region of 130 DQV lines, a 
peripheral annulus of 140 SQV lines, and the exterior 
vortex-free CF region of equilibrium width. The bound- 
ary between the DQV and SQV regions is relatively sharp 
(10 % of the radius of the array, or additionally roughly 
40 DQV and 50 SQV lines). 

The coaxial configuration is stable against externally 
applied perturbations. In fact, no means were found 
to change the configuration, without removing vortex 
lines by deceleration or by adding new ones during ac- 
celeration. The noise in the rotation drive is of order 
|Af2| ~ 1 • 10~ 3 rad/s. However, no noticeable changes 
in the coaxial configuration were observed when the ar- 
ray is maintained for several hours at constant f2. This 
is not surprising since even a sinusoidal modulation of Q 
with an amplitude Af2 ~ 0.1 rad/s does not cause mea- 
surable changes. (The modulation by 2Af2 is less than 
the separation between the annihilation and critical ve- 
locity thresholds and thus the number of vortex lines is 
not yet changed.) This suggests that any spontaneous re- 
ordering of the two vortex species within the array from 
the segregated domain structure to a mixed configuration 
is below our measuring resolution. Thus the cool down 
from T c to the measuring temperature or a magnetic field 
sweep after the transition through T c are not expected to 
introduce changes in the array configuration. 

The segregation into two coaxial domains, with the 
SQV lines in the outer annulus, is always observed when 
both vortex species are found to coexist after cooling 
through T c . We also find, by measuring the width of 
the CF annulus |l7) , that the total number of lines ad- 
justs itself so that the energy in Eq. (^j|) is minimized. 
This is expected, since the total circulation corresponds 
to a dominant energy term. In comparison the structure 
of an individual vortex line or the spatial distribution of 
the two vortex species within the array, represent much 
smaller energy contributions. 

C. Sequentially formed vortex cluster 

The stability of vortex arrays in superfluid 3 He can also 
be approached from the other extreme, by studying the 
disorder in artificially constructed arrays of coexisting 
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FIG. 14. Radial distributions N±(r 2 ) and N 2 (r 2 ) in a 
sequentially formed vortex cluster. The original array was 
an equilibrium state of SQV lines, to which DQV lines were 
added later by increasing rotation in the superfluid state. The 
SQV array was prepared at 29.3 bar pressure and 9.91 mT 
field, by cooling through T c at 7 /iK/min and 0.6 rad/s. Cool- 
ing was continued until 0.81 T c , where Q. was slowly increased 
to 1.4 rad/s, to create the DQV lines. The final step was a 
slow deceleration to zero, during which the radial distribu- 
tions Ni(r 2 ) and N2(r 2 ) were recorded. Both Q sweeps were 
performed at \dQ/dt\ = 2 • 10~ 3 rad/s" 2 . The critical CF ve- 
locity of vortex formation at 1.4 rad/s was £l c — 0.24 rad/s, 
which corresponds to a linear velocity of v c = 0.60 mm/s. 
The inset illustrates schematically the array configuration, a 
plastically deformed agglomerate of vortex lines. 



SQV and DQV lines. These can be prepared by making 
use of the metastability of existing vorticity. 

We use the following procedure: The sample is first 
cooled through T c at £1 < 0.6 rad/s and H > H c , so 
that an equilibrium vortex state with only SQV lines 
is formed. DQV lines are added later, by increasing 
slowly while the temperature is already below T c , such 
that the desired number of DQV lines is created at a 
stable critical counterflow velocity v c = Q C R Hj. In the 
presence of an existing array this value of CF represents 
the separation between the annihilation and critical ve- 
locity thresholds. The DQV lines are here injected one 
at a time into the existing array, presumably from the 
same nucleation center on the cylindrical wall where the 
superflow instability takes place (i.e. they are thought to 
emanate from a rectilinear source which is close to the 
cylindrical wall and parallel to the rotation axis 0). 

In Fig. [li] the radial decomposition of an array is shown 
which has been prepared with this procedure, by adding 
150 DQV lines to an existing equilibrium array of 310 
SQV lines. During deceleration, initially both types of 
vortices are seen to annihilate, but on moving closer to 
the center the proportion of DQV lines is reduced and the 
SQV lines start to dominate. Such an annihilation record 
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could correspond to the configuration shown schemati- 
cally in the inset of Fig. [lj, which in turn could result 
from the process by which the array was formed. The 
DQV lines are here formed in a regular periodic process, 
where a newly created DQV line arrives roughly to the 
same spot at the edge of the array. It makes space for 
itself by pushing the existing lines aside. In this way the 
DQV lines end up residing in a segregated pocket at the 
edge of the cluster, if no mixing of lines occurs later. 

From the data in Fig. |l4] it is not possible to judge 
how sharp the boundary is between the two domains and 
to what extent mixing of the two vortex species might 
have occurred. However, no significant change from this 
configuration is noticed even if the cluster is subjected to 
rapid modulation of the rotation velocity with an am- 
plitude which does not yet add new vortices or anni- 
hilate existing ones. On this basis we expect sponta- 
neous reordering or annealing not to be of importance. 
The sequentially formed array is most likely a state with 
much reduced long-range order, which resembles a ran- 
dom glassy agglomerate, but is subject to hexagonal 
nearest-neighbor coordination and a circular boundary 
condition. 



D. Homogeneity of superfluid transition 

In Chap. we have focused on two questions: (1) 
the metastability and width of the SQV <-* DQV tran- 
sition and (2) the configuration of the vortex array in 
the coexistence regime. We might now simply conclude 
that the segregated coaxial array, which is observed af- 
ter a cool down through T c in the coexistence regimes, is 
the expected minimum energy configuration. Since vor- 
tex clusters in 3 He-A cannot relax during the course of 
the experiments, this configuration must also be the one 
which was originally formed in the transition through T c . 

However, in any real experiment the superfluid transi- 
tion will acquire some degree of inhomogeneity. On some 
level this will lead to irreversibility and will interfere with 
the argument that the transition is adiabatic and mini- 
mizes the total energy. Here the main source for inho- 
mogeneity is the thermal resistance concentrated in the 
orifice of the container. The liquid 3 He volume consists 
of the sample container and the sintered heat exchanger 
in contact with the refrigerator. These are connected to- 
gether with a long tube, which ends in the orifice on the 
bottom of the sample cylinder |18|. During cool down a 
thermal gradient is present, in addition to residual heat 
leaks, which maintain the top end of the sample container 
as the warmest point in the 3 He column. This means that 
the heat flow through the orifice is composed of the cool- 
ing rate of the liquid and the residual heat leak. The 
residual heating consists of the NMR absorption and of 
thermal relaxation losses seeping out of the walls of the 
sample container. 



At a cooling rate of 1 /iK/min the heat flow from the 
container is of order 0.5 nW and the thermal gradient in- 
side the container ~ 0.5 /iK/mm (at T c ~ 2.4 mK). The 
superfluid transition emerges thus as a transition front 
through the orifice on the bottom of the cylinder. Ini- 
tially it adopts a hemispherical shape B. The width of 
the phase front is determined by the local thermal gra- 
dient and thermal fluctuations. Its velocity can be con- 
trolled externally by adjusting the cooling rate. The heat 
flow, in turn, is mainly carried by a superfluid thermal 
countercurrent, which appears below T c and converges 
towards the orifice. It is oriented predominantly in the 
axial direction, transverse to the rotational flow. 

Clearly this picture of the transition is not that of a ho- 
mogeneous second order transition (although it is proba- 
bly closer to it than in any other macroscopic system) . In 
equilibrium conditions at low the SQV is the minimum 
energy vortex structure, but thermal CF or other inho- 
mogeneities provide a bias for DQV lines to be formed. 
This perturbation is concentrated in the neighborhood of 
the orifice, from where vortex formation starts during the 
superfluid transition. Here DQV lines are created first 
and SQV lines appear later when the transition front ex- 
pands into the more homogeneous peripheral region of 
the cylinder. This interpretation, which is suggested by 
the dependence of the coexistence regime on the cool- 
ing rate at T c (Fig. |l^) , provides a second reason for 
the coaxial array configuration with the SQV lines in the 
outer annulus. 

According to this interpretation, the interface between 
the two coaxial domains becomes a similar phase bound- 
ary as that in Fig. O, but now as a function of the radial 
coordinate r of the container, along which the cooling 
rate varies as the normal-to-superfluid phase front moves 
through the cylinder. This alternative explanation thus 
rests on the assumption that the cooling rate at T c is one 
of the maj or controlling factors of the phase boundary 
in Fig. [12|: SQV lines are formed closer to equilibrium 
conditions while nonequilibrium and perturbations favor 
the formation of DQV lines. 



VI. CONCLUSION 

By comparing experiment to numerical simulation, we 
conclude that below T c in superfluid 3 He vortex arrays 
are formed by plastic deformation: Vortex lines are frozen 
in the array to the configuration in which the array was 
formed. High energy barriers separate different array 
configurations and thermally activated processes (or any 
other mechanisms) are not able to anneal faults, to cre- 
ate a new structure with lower energy and improved long 
range order. 

In arrays formed by cooling through T c in rotation, 
a segregated coaxial configuration is observed, with the 
SQV lines located in the outer annulus. Our numerical 
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analysis of 2D arrays with discrete rectilinear vortex lines 
confirms that this coaxial configuration is energetically 
the most preferable state in the coexistence regime. It 
is distinguished from all other segregated domain struc- 
tures by a gain in the outer surface energy. This energy 
reduction arises on the scale of the intervortex distance 
rather than the container dimensions, due to a smoother 
outer interface, when the border is formed by vortex lines 
with a smaller quantum number. 
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